Pilot-scale genome-wide association mapping in diverse sorghum germplasms identified novel genetic loci linked to major agronomic, root and stomatal traits.

This genome-wide association studies (GWAS) used a subset of 96 diverse sorghum accessions, constructed from a large collection of 219 accessions for mining novel genetic loci linked to major agronomic, root morphological and physiological traits. The subset yielded 43,452 high quality single nucleotide polymorphic (SNP) markers exhibiting high allelic diversity. Population stratification showed distinct separation between caudatum and durra races. Linkage disequilibrium (LD) decay was rapidly declining with increasing physical distance across all chromosomes. The initial 50% LD decay was ~ 5 Kb and background level was within ~ 80 Kb. This study detected 42 significant quantitative trait nucleotide (QTNs) for different traits evaluated using FarmCPU, SUPER and 3VmrMLM which were in proximity with candidate genes related and were co-localized in already reported quantitative trait loci (QTL) and phenotypic variance (R2) of these QTNs ranged from 3 to 20%. Haplotype validation of the candidate genes from this study resulted nine genes showing significant phenotypic difference between different haplotypes. Three novel candidate genes associated with agronomic traits were validated including Sobic.001G499000, a potassium channel tetramerization domain protein for plant height, Sobic.010G186600, a nucleoporin-related gene for dry biomass, and Sobic.002G022600 encoding AP2-like ethylene-responsive transcription factor for plant yield. Several other candidate genes were validated and associated with different root and physiological traits including Sobic.005G104100, peroxidase 13-related gene with root length, Sobic.010G043300, homologous to Traes_5BL_8D494D60C, encoding inhibitor of apoptosis with iWUE, and Sobic.010G125500, encoding zinc finger, C3HC4 type domain with Abaxial stomatal density. In this study, 3VmrMLM was more powerful than FarmCPU and SUPER for detecting QTNs and having more breeding value indicating its reliable output for validation. This study justified that the constructed subset of diverse sorghums can be used as a panel for mapping other key traits to accelerate molecular breeding in sorghum.

Parallel tuning of semi-dwarfism via differential splicing of Brachytic1 in commercial maize and smallholder sorghum.

In the current genomic era, the search and deployment of new semi-dwarf alleles have continued to develop better plant types in all cereals. We characterized an agronomically optimal semi-dwarf mutation in Zea mays L. and a parallel polymorphism in Sorghum bicolor L. We cloned the maize brachytic1 (br1-Mu) allele by a modified PCR-based Sequence Amplified Insertion Flanking Fragment (SAIFF) approach. Histology and RNA-Seq elucidated the mechanism of semi-dwarfism. GWAS linked a sorghum plant height QTL with the Br1 homolog by resequencing a West African sorghum landraces panel. The semi-dwarf br1-Mu allele encodes an MYB transcription factor78 that positively regulates stalk cell elongation by interacting with the polar auxin pathway. Semi-dwarfism is due to differential splicing and low functional Br1 wild-type transcript expression. The sorghum ortholog, SbBr1, co-segregates with the major plant height QTL qHT7.1 and is alternatively spliced. The high frequency of the Sbbr1 allele in African landraces suggests that African smallholder farmers used the semi-dwarf allele to improve plant height in sorghum long before efforts to introduce Green Revolution-style varieties in the 1960s. Surprisingly, variants for differential splicing of Brachytic1 were found in both commercial maize and smallholder sorghum, suggesting parallel tuning of plant architecture across these systems.

Dissecting Adaptive Traits with Nested Association Mapping: Genetic Architecture of Inflorescence Morphology in Sorghum.

In the cereal crop sorghum (Sorghum bicolor) inflorescence morphology variation underlies yield variation and confers adaptation across precipitation gradients, but its genetic basis is poorly understood. We characterized the genetic architecture of sorghum inflorescence morphology using a global nested association mapping (NAM) population (2200 recombinant inbred lines) and 198,000 phenotypic observations from multi-environment trials for four inflorescence morphology traits (upper branch length, lower branch length, rachis length, and rachis diameter). Trait correlations suggest that lower and upper branch length are under somewhat independent control, while lower branch length and rachis diameter are highly pleiotropic. Joint linkage and genome-wide association mapping revealed an oligogenic architecture with 1-22 QTL per trait, each explaining 0.1-5.0% of variation across the entire NAM population. There is a significant enrichment (2.twofold) of QTL colocalizing with grass inflorescence gene homologs, notably with orthologs of maize Ramosa2 and rice Aberrant Panicle Organization1 and TAWAWA1 Still, many QTL do not colocalize with inflorescence gene homologs. In global georeferenced germplasm, allelic variation at the major inflorescence QTL is geographically patterned but only weakly associated with the gradient of annual precipitation. Comparison of NAM with diversity panel association suggests that naive association models may capture some true associations not identified by mixed linear models. Overall, the findings suggest that global inflorescence diversity in sorghum is largely controlled by oligogenic, epistatic, and pleiotropic variation in ancestral regulatory networks. The findings also provide a basis for genomics-enabled breeding of locally-adapted inflorescence morphology.

Interaction Between Induced and Natural Variation at oil yellow1 Delays Reproductive Maturity in Maize.

We previously demonstrated that maize (Zea mays) locus very oil yellow1 (vey1) encodes a putative cis-regulatory expression polymorphism at the magnesium chelatase subunit I gene (aka oil yellow1) that strongly modifies the chlorophyll content of the semi-dominant Oy1-N1989 mutants. The vey1 allele of Mo17 inbred line reduces chlorophyll content in the mutants leading to reduced photosynthetic output. Oy1-N1989 mutants in B73 reached reproductive maturity four days later than wild-type siblings. Enhancement of Oy1-N1989 by the Mo17 allele at the vey1 QTL delayed maturity further, resulting in detection of a flowering time QTL in two bi-parental mapping populations crossed to Oy1-N1989 The near isogenic lines of B73 harboring the vey1 allele from Mo17 delayed flowering of Oy1-N1989 mutants by twelve days. Just as previously observed for chlorophyll content, vey1 had no effect on reproductive maturity in the absence of the Oy1-N1989 allele. Loss of chlorophyll biosynthesis in Oy1-N1989 mutants and enhancement by vey1 reduced CO2 assimilation. We attempted to separate the effects of photosynthesis on the induction of flowering from a possible impact of chlorophyll metabolites and retrograde signaling by manually reducing leaf area. Removal of leaves, independent of the Oy1-N1989 mutant, delayed flowering but surprisingly reduced chlorophyll contents of emerging leaves. Thus, defoliation did not completely separate the identity of the signal(s) that regulates flowering time from changes in chlorophyll content in the foliage. These findings illustrate the necessity to explore the linkage between metabolism and the mechanisms that connect it to flowering time regulation.

Genetic Architecture of Chilling Tolerance in Sorghum Dissected with a Nested Association Mapping Population.

Dissecting the genetic architecture of stress tolerance in crops is critical to understand and improve adaptation. In temperate climates, early planting of chilling-tolerant varieties could provide longer growing seasons and drought escape, but chilling tolerance (<15°) is generally lacking in tropical-origin crops. Here we developed a nested association mapping (NAM) population to dissect the genetic architecture of early-season chilling tolerance in the tropical-origin cereal sorghum (Sorghum bicolor [L.] Moench). The NAM resource, developed from reference line BTx623 and three chilling-tolerant Chinese lines, is comprised of 771 recombinant inbred lines genotyped by sequencing at 43,320 single nucleotide polymorphisms. We phenotyped the NAM population for emergence, seedling vigor, and agronomic traits (>75,000 data points from ∼16,000 plots) in multi-environment field trials in Kansas under natural chilling stress (sown 30-45 days early) and normal growing conditions. Joint linkage mapping with early-planted field phenotypes revealed an oligogenic architecture, with 5-10 chilling tolerance loci explaining 20-41% of variation. Surprisingly, several of the major chilling tolerance loci co-localize precisely with the classical grain tannin (Tan1 and Tan2) and dwarfing genes (Dw1 and Dw3) that were under strong directional selection in the US during the 20th century. These findings suggest that chilling sensitivity was inadvertently selected due to coinheritance with desired nontannin and dwarfing alleles. The characterization of genetic architecture with NAM reveals why past chilling tolerance breeding was stymied and provides a path for genomics-enabled breeding of chilling tolerance.

Adult plant resistance in maize to northern leaf spot is a feature of partial loss-of-function alleles of Hm1.

Adult plant resistance (APR) is an enigmatic phenomenon in which resistance genes are ineffective in protecting seedlings from disease but confer robust resistance at maturity. Maize has multiple cases in which genes confer APR to northern leaf spot, a lethal disease caused by Cochliobolus carbonum race 1 (CCR1). The first identified case of APR in maize is encoded by a hypomorphic allele, Hm1A, at the hm1 locus. In contrast, wild-type alleles of hm1 provide complete protection at all developmental stages and in every part of the maize plant. Hm1 encodes an NADPH-dependent reductase, which inactivates HC-toxin, a key virulence effector of CCR1. Cloning and characterization of Hm1A ruled out differential transcription or translation for its APR phenotype and identified an amino acid substitution that reduced HC-toxin reductase (HCTR) activity. The possibility of a causal relationship between the weak nature of Hm1A and its APR phenotype was confirmed by the generation of two new APR alleles of Hm1 by mutagenesis. The HCTRs encoded by these new APR alleles had undergone relatively conservative missense changes that partially reduced their enzymatic activity similar to HM1A. No difference in accumulation of HCTR was observed between adult and juvenile plants, suggesting that the susceptibility of seedlings derives from a greater need for HCTR activity, not reduced accumulation of the gene product. Conditions and treatments that altered the photosynthetic output of the host had a dramatic effect on resistance imparted by the APR alleles, demonstrating a link between the energetic or metabolic status of the host and disease resistance affected by HC-toxin catabolism by the APR alleles of HCTR.

Population genomics of sorghum (Sorghum bicolor) across diverse agroclimatic zones of Niger.

Improving adaptation of staple crops in developing countries is important to ensure food security. In the West African country of Niger, the staple crop sorghum (Sorghum bicolor) is cultivated across diverse agroclimatic zones, but the genetic basis of local adaptation has not been described. The objectives of this study were to characterize the genomic diversity of sorghum from Niger and to identify genomic regions conferring local adaptation to agroclimatic zones and farmer preferences. We analyzed 516 Nigerien accessions for which local variety name, botanical race, and geographic origin were known. We discovered 144 299 single nucleotide polymorphisms (SNPs) using genotyping-by-sequencing (GBS). We performed discriminant analysis of principal components (DAPC), which identified six genetic groups, and performed a genome scan for loci with high discriminant loadings. The highest discriminant coefficients were on chromosome 9, near the putative ortholog of maize flowering time adaptation gene Vgt1. Next, we characterized differentiation among local varieties and used a genome scan of pairwise FST values to identify SNPs associated with specific local varieties. Comparison of varieties named for light- versus dark-grain identified differentiation near Tannin1, the major gene responsible for grain tannins. These findings could facilitate genomics-assisted breeding of locally adapted and farmer-preferred sorghum varieties for Niger.

Increased Power To Dissect Adaptive Traits in Global Sorghum Diversity Using a Nested Association Mapping Population.

Adaptation of domesticated species to diverse agroclimatic regions has led to abundant trait diversity. However, the resulting population structure and genetic heterogeneity confounds association mapping of adaptive traits. To address this challenge in sorghum [Sorghum bicolor (L.) Moench]-a widely adapted cereal crop-we developed a nested association mapping (NAM) population using 10 diverse global lines crossed with an elite reference line RTx430. We characterized the population of 2214 recombinant inbred lines at 90,000 SNPs using genotyping-by-sequencing. The population captures ∼70% of known global SNP variation in sorghum, and 57,411 recombination events. Notably, recombination events were four- to fivefold enriched in coding sequences and 5' untranslated regions of genes. To test the power of the NAM population for trait dissection, we conducted joint linkage mapping for two major adaptive traits, flowering time and plant height. We precisely mapped several known genes for these two traits, and identified several additional QTL. Considering all SNPs simultaneously, genetic variation accounted for 65% of flowering time variance and 75% of plant height variance. Further, we directly compared NAM to genome-wide association mapping (using panels of the same size) and found that flowering time and plant height QTL were more consistently identified with the NAM population. Finally, for simulated QTL under strong selection in diversity panels, the power of QTL detection was up to three times greater for NAM vs. association mapping with a diverse panel. These findings validate the NAM resource for trait mapping in sorghum, and demonstrate the value of NAM for dissection of adaptive traits.

Comparative Transcriptome and Lipidome Analyses Reveal Molecular Chilling Responses in Chilling-Tolerant Sorghums.

Chilling temperatures (0 to 15°C) are a major constraint for temperate cultivation of tropical-origin crops, including the cereal crop sorghum ( [L.] Moench). Northern Chinese sorghums have adapted to early-season chilling, but molecular mechanisms of chilling tolerance are unknown. We used RNA sequencing of seedlings to compare the chilling-responsive transcriptomes of a chilling-tolerant Chinese accession with a chilling-sensitive US reference line, and mass spectrometry to compare chilling-responsive lipidomes of four chilling-tolerant Chinese accessions with two US reference lines. Comparative transcriptomics revealed chilling-induced up-regulation of cold-response regulator C-repeat binding factor (CBF) transcription factor and genes involved in reactive oxygen detoxification, jasmonic acid (JA) biosynthesis, and lipid remodeling phospholipase Dα1 (α) gene in the chilling-tolerant Chinese line. Lipidomics revealed conserved chilling-induced increases in lipid unsaturation, as well as lipid remodeling of photosynthetic membranes that is specific to chilling-tolerant Chinese accessions. Our results point to CBF-mediated transcriptional regulation, galactolipid and phospholipid remodeling, and JA as potential molecular mechanisms underlying chilling adaptation in Chinese sorghums. These molecular systems underlying chilling response could be targeted in molecular breeding for chilling tolerance.